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PHASE VARIATION IN THE PULSE PROFILE OF SMC X-l 
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ABSTRACT 

We present the results of the timing and spectral analysis of X-ray high-state observations of the high-mass 
X-ray pulsar SMC X-l with Chandra, XMM-Newton, and ROSAT, taken between 1991 and 2001. The source 
hasLx ~ (3-5) x 10 38 ergs s" 1 , and the spectra can be modeled as a power law plus blackbody with kT^B ~ 0.18 
keV and a reprocessed emission radius Rbb ~2x 10 s cm, assuming a distance of 60 kpc to the source. Energy- 
resolved pulse profiles show several distinct forms, more than half of which include a second pulse in the soft 
profile, previously documented only in hard energies. We also detect significant variation in the phase shift 
between hard and soft pulses, as has recently been reported in Her X-l. We suggest an explanation for the 
observed characteristics of the soft pulses in terms of precession of the accretion disk. 

Subject headings: accretion, accretion disks — stars: neutron — X-rays: binaries: individual (SMC X-l) — 
stars: pulsars: individual (SMC X-l) 

1. INTRODUCTION 

SMC X-l is a high-mass X-ray binary system consisting 
of an X-ray source, a 1 . 6 pulsar, accret ing from a 17.2 
Mq companion, Sk 160 (Clarkson et al. 2003). It is also the 
only X-ray puls ar for which no spin-down episodes have been 
observed (Kahabka & Li 1999). The X-ray pulsar has a spin 
period of 0.71 s and an eclipsing orbital period of 3.89 days, 
and in addition it shows a superorbital variation in flux with 
period between 45 and 60 days. This is likely due to p reces- 
sion of a warped accretion disk ( Woidowski et all!998l) . 

The presence of the accretion disk is likely to have signif- 
icant consequences for the observed pulse profile. In Her 
X-l, a low-mass X-ray pulsar with a regular 35 day super- 
orbital period, reprocessing of hard X-rays by the disk give s 
rise to a pu l sating , soft spectral component (Endo et al. 2000). 
IZane et aTJ J2004I) have shown that the phase difference be- 
tween hard and soft pulses changes as the accretion disk pre- 
cesses. If simil ar reprocessin g is at work in SMC X-l, as ar- 
gued by Hickox et alj J2004I) . the pulse profiles should show 
a variation similar to that of Her X-l. To that end, we have 
analyzed pulse profiles of SMC X-l in different epochs from 
Chandra, XMM, and ROSAT. In § 2 we discuss the observa- 
tions and our methods of data analysis. In § 3 we present the 
results of our analysis, and in § 4 we consider the implications 
of our results and future work. 

2. OBSERVATIONS 

The observations reported in this Letter are described in 
Table 1 and were made with Chandra's ACIS-S in continu- 
ous clocking (CC) mode (CXC Proposers' Observatory Guide 
v5.0, 2002), the EPIC-pn detectors on XMM-Newton in full 
window (X101) and small window (X201) modes (XMM- 
Newton Users' Handbook v2.2, 2004), and the ROSAT PSPC 
detectors in pointing mode (ROSAT Users' Handbook 1996). 
We shall refer to the observations as C104, X201, RnOO, etc. 
(see Table 1). Fig. 1 locates the Chandra and XMM observa- 
tions in the All-Sky Monitor (ASM) lightcurve from the Rossi 
X-ray Timing Explorer (RXTE) mission. We estimate the su- 

1 Kenyon College, Department of Physics, Gambier, OH 43022; 
neilsenj@kenyon.edu 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cam- 
bridge, MA 02138; rhickox@cfa.harvard.edu; saku@head.cfa.harvard.edu 



FIG. 1. — RXTE ASM lightcurve for SMC X-l with arrows designating sev- 
eral observations. Quick-look results as provided by the RXTE/ASM team. 



perorbital phase (</>so) as 0.16 for the Chandra observations, 
0.37 for X101, and 0.42 for X201, where S o = at the be- 
ginning of the X-ray high state. Given the variation in (f>so, 
the uncertainty for these estimates is plus or minus a few days 
or 4>so ~ 0.05. The 0so estimates for ROSAT are unavailable 
because those observations were taken before the launch of 
RXTE. 

Since the CC mode provides only one dimension of spa- 
tial information, ACIS source events were extracted from 
a source-centered rectangle of size 0."492 x 1 ."968. Back- 
ground events were extracted from two regions 0."492 x 9. "84 
equidistant from the source. For X101, which shows the 
source entering an eclipse, events were located within a cir- 
cular region of radius 45" and extracted from the 4.9 ks of 
pre-eclipse data; for X201, we used a circular region of radius 
51 "2. We ignored the XMM background because its surface 
brightness was less than 0.1% of that of the source region. 
ROSAT source events were extracted from a circular region of 
radius 2. '5. For Ra02, which shows an eclipse, events were 
taken from the 5.5 ks of post-eclipse data. 

3. ANALYSIS 

Before the spectroscopic and timing analysis, we performed 
a barycentric correction and corrected for the orbital mo- 
tion of SMC X-l using the orbital ephemeris calculated by 
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TABLE 1 
Observations of SMC X- 1 



Observation 


Observation 


Mission 


Start Time 


Observation 


Period 


ID Number 


Reference 




(JD-2400000) 


Time (ks) 


(s) 


rp400022nOO 


RnOO 


ROSAT 


48536.67 


16.6 


0.709114(2) 


rp400022a02 


Ra02 


ROSAT 


49141.50 


11.8 


0.708598(5) 


400102 


C102 


Chandra 


51832.17 


6.2 


0.70567(2) 


400103 


C103 


Chandra 


51833.34 


6.1 


0.70567(2) 


400104 


C104 


Chandra 


51834.31 


6.5 


0.70567(1) 


0011450101 


X101 


XMM 


52060.59 


46.4 


0.70542(2) 


0011450201 


X201 


XMM 


52229.65 


40.0 


0.70522(3) 



Woidowsk i et al.l (^998). The following analysis was per- 
formed with CIAO version 3.1 for Chandra and ROSAT data 
and SAS version 5.3.3 for XMM. 

3.1. Phase-Averaged Spectroscopy 

For Chandra and XMM, we extracted phase-averaged spec- 
tra, modeling the 0.6-9.8 keV emissions with neutral absorp- 
tion, a blackbody with kTss ~ 0.18 keV, and a power law. 
For the XMM spectra we also included a high- energy cutoff 
(E cut - 6 keV, Eton ~ 8 keV lWoo etali n^). For Chan- 
dra and X201, we excluded the energy range 1.7-2.8 keV to 
avoid instrumental effects (near the Au and Si edges) for high 
count data (Mil ler et alJ l2002). Although the excellent time 
resolution of CC mode prevents some photon pileup, w e used 
the IS IS pileup kernel for the Chandra observations JPavisI 
1MB, which show a pileup fraction of ^17%. For XMM we 
minimized pileup in the spectra by excluding events from the 
central one-third radius of the extraction regions. 

Spectral fits from C104, X201, and X101 are shown in Fig. 
2 and results are listed in Table 2. All spectra are dominated 
by blackbody emission below ~1.0 keV, athough the similar 
wavy residuals in the Chandra and X201 spectra below 1.0 
keV suggest that the soft emission is not exactly blackbody. 
For X101, the data are taken as the source is entering eclipse, 
and the spectrum could not be fitted using the simple model 
above. Fixing T to a typical value of 0.9, we found that a 
good fit is achieved by including a partial covering absorp- 
tion, with a covering fraction of 88%. This accounts for the 
extra absorption by the dense gas around the star as the eclipse 
begins, but allows for some X-rays to be scattered around the 
absorbing region. Residuals suggest that there may also be Fe 
absorption at ^6.5 keV, but more detailed spectral analysis 
is beyond the scope of this Letter. Since blackbody emission 
dominates below 1.0 keV, and power-law emission above 2.0 
keV, we can use energy-resolved profiles in these ranges to 
examine the variation of the separate components. More de- 
tailed analysis, using pulse-phase spectroscopy, will be pre- 
sented in a forthcoming paper. 

3.2. Timing Analysis 

Since SMC X-l is a fast-period pulsar, pulse profile analy- 
sis requires a time resolution of better than 0.7 1 s. Fortunately, 
the Chandra, X201, X101, and ROSAT observations have suf- 
ficient time resolutions of 2.85, 6, 73.4, and 130 ms, respec- 
tively. We used ef search in the software package XRONOS to 
measure the pulse periods (see Tabl e 1). These are consistent 
with the spin-up rate calculate d bv IWoi dowski et alJ J1998I) 
and the periods determined by Vrtilek et al. (2001) from the 
same Chandra data. We extracted high- and low-energy pulse 
profiles from each of the seven observations (see Fig. 3). Here 




1 2 

channel energy (keV) 




channel energy (keV) 

FIG. 2.— Spectra from (a) C104, (b) X201, and (c) X101 with blackbody 
and power law components shown. As X101 is a pre-eclipse observation, it 
was fit with an additional partial covering absorption. Note that in (a) the 
model components do not add exactly to the observed counts. This is due to 
the pileup kernel used in the fit. 



PULSE PROFILE OF SMC X-l 
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TABLE 2 
Fits to SMC X-l Spectra 



l cimm c ter 


CI 02 


CI 03 


CI 04 


X201 


X101 




°- 14 toi 

0.89*| 
n 177+0 007 

ua/ '-0.007 


23+ 02 

CP 




°- 17+ oo! 




r 

kT BB 


1 03+M3 
-0 04 

169*™ 5 

u - 1Dy -0.005 


0.92*| 
n no+0.005 
U - ,/S -0.005 


0.91*1 
0.190|p 


0.9 
173+ 0013 

8 9 +IT 


Beat (keV) 
Efoid (keV) 










Pileup fraction 


0.16 


0.18 


0.16 






Partial N H (10 22 cirr 2 ) 










1 2+ ' 1 


Covering fraction 










0.88+^ 2 


XI (d.o.f.) 


1.36 (384) 


1.19(387) 


1.22 (390) 


1.20(516) 


1.20 (258) 


Observed flux (10 9 ergs cm 2 s ) 


1.1 


1.1 


1.1 


1.0 


0.57 


Unabsorbed flux (1CT 9 ergs cm -2 s _1 ) 


1.2 


1.2 


1.3 


1.1 


0.84 


Luminosity (ergs s ) 


5.2 


5.3 


5.4 


4.7 


3.6 



NOTE. — Errors are 90% confidence for a single parameter. Fluxes are for 0.6-9.8 keV. A distance of 60 kpc is assumed for the 
SMC. 



Normalized counts/s 




FIG. 3. — Hard and soft pulse profiles for (a) ROSAT, (b) Chandra, and (c) 
XMM. The hard range is 2.0-8.0 keV (1.5-2.4 keV for ROSAT); soft is 0.5- 
1 .0 keV. Two pulse periods are shown for clarity. Note that the pulse phases 
shown are not absolute, but are only relative to each observation. 



we do not correct for pileup effects, since these have minimal 
impact on the shapes of the pulses. 

All seven observations show the double-peaked hard 
pulse profile th at has been thor oughly documented 
(Woidowsk ietafl fl998t iPaul et alJ 12003 iNaik & Paull 
120041) . The soft profile, however, varies markedly. While 
earlier studies of S MC X-l have shown a roughly sinu soidal 
shape of soft pulses (IPaul et aU2002t INaik & Paull2004l) . the 
ROSAT, Chandra, and XMM observations show a complex 
variation of profiles. Although the ROSAT pulse profiles 
show a single main soft peak, four of the five Chandra and 
XMM observations show double-peaked soft pulses. 

The first ROSAT observation, RnOO, shows a single asym- 
metric hump at lower energies, slightly out of phase with the 
hard pulses. For Ra02, the soft pulses are lightly double- 
peaked, although the poor spectral resolution of the PSPC 
places the exact details of the profiles in doubt. The X101 
pulse profile shows a single broad peak quite similar to RnOO. 
In the Chandra data, we see a clear second peak in all three 
soft pulse profiles, and we note a significant phase difference 
between hard and soft pulses. In X201 we also find double- 
peaked soft pulses, but these are almost in phase with the hard 
pulses. 

For Chandra and XMM we cross-correlated the hard and 
soft profiles, to quantify the phase sh ift between the pulses 
(see Fig. 4), as in Ramsav et al. (2002) and Zane et al. (2004). 
The Chandra profiles show strong anticorrelations at a phase 
shift of 0° and positive correlations at ±90°. We see mul- 
tiple peaks because the hard and soft profiles each have two 
pulses. Observation X101 shows a broad correlation with a 
peak at ^20°; in X201 we detect possible phase shifts of ^0°, 
-±90°, and -±170°. 

4. DISCUSSION 

In these observations both hard and soft X-rays show pul- 
sations, but differences in the hard and soft pulse profiles in- 
dicate a dif ferent geometric or physical origin of emission 
iPaul et al] 2002). Soft pulses in the most luminous X-ray pul- 
sars are likely to originate at the inner edge of the accretion 
disk, where hard pulses from t he neutron star are reradiated at 
lower energies by disk gases (Hickox et al. 2004). If this pic- 
ture is valid for SMC X-l, the observed blackbody component 
should be consistent with emission from the inner disk. For 
reprocessing by a partial spherical surface around the neutron 
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FIG. 4. — Cross-correlation between hard and soft pulse profiles for Chan- 
dra and XMM data. 



star, the blackbody radius Rbb is given bv JPaul et all2002l) 



Rrr — 



Lx 



BB 



4 x 10 38 ergs s we have 



Using kT BB ~ 0.18 keV and L x 
fl B B ~ 1.7xl0 8 cm. 

In the standard model for X-ray pulsars, the inner disk ra- 
dius /?o is close to the magnetospheric radius R,„, where the 
magnetic pressure of the dipole field equals the ram pressure 
of the infalling gas. However, R m is difficult to estimate for 
SMC X-l, because the surface fi-field has not been directly 
measured. We note that R B b is close to the corotation ra- 
dius R cm = (GMP 2 /4ir 2 ) 1 / 3 ~ 1.3 x 10 8 cm, which is another 
estimate for Ro. However, the lack o f spin-down ep isodes 
for SMC X-l suggests that R < R col - (Kahabka &l1|1999I) . 
The situation is therefore unclear, but without additional con- 
straints it is certainly possible that the soft pulses originate at 
the inner accretion disk. 

In light of this picture we consider the following pulse pro- 
file characteristics: 

1. The hard pulse profile is consistently double-peaked, 
typically with one dominant peak. 

2. The soft pulse profile varies in shape, and has been 
observed with either a single asymmetric peak or two 
peaks with varying degrees of symmetry. 



The pulse profiles exhibit large-scale variations in 
phase relationship. Closely-spaced observations have 
revealed identical phase shifts. 



IZane et all fe00 4) suggest that the varying hard-soft phase 
shift in Her X-l may be due to the precession of the accre- 
tion disk, so we consider the possibility that such a disk could 
reproduce the above characteristics in SMC X-l. If the ob- 
servational line of sight is sufficiently close to the plane of the 
magnetic axis, we might always see two hard pulses, one from 
each magnetic pole. However, the soft profile could change 
dramatically in shape and relative phase, as observed, if pre- 
cession of the disk causes the visible part of the reprocessing 
region to vary. 

We note however that we cannot conclusively verify this 
picture and its analogy to Her X-l. First, we lack sufficient 
sampling of the full 60 day superorbital period. Second, the 
hard-soft phase shift appears to vary more rapidly for SMC 
X-l (a change of ~90° between SO ~ 0.16 and 0.42) than 
for H er X-l (-20° between S o ^ 0.03 and 0.17 lZane et ail 
120041) . Third, it is not completely clear even in Her X-l that 
the phase variation is caused by the 35 day disk precession; 
soon after the observations of Zane et al. (2004), the source 
entered an anomalous low state, an event that may be caused 
by changes in the inne r disk structure JBovd & Stilll 120041 
Oosterbroek et al 1l200Th . Such changes could alter the shape 
of the inner reprocessing region, and thus the soft profiles, in 
a way that is not directly related to the superorbital period. 

We conclude that in SMC X-l, we have observed variations 
in the soft pulse profile that likely reflect changes in the repro- 
cessing of hard X-rays by the accreting gas. These variations 
may be related to the superorbital precession of the accretion 
disk, as has been proposed for Her X- 1 . If it can be shown 
with more observations that the pulse profiles and superor- 
bital period are correlated, then we will have the opportunity 
to apply significant constraints to the geometry of SMC X-l. 
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